Landauer-like formula for dissipative tunneling 
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The Landauer formula for electrical conductance is simple but works remarkably well in meso- 
scopic systems. We propose a Landauer-like formula for calculating an escape rate out of a dissipative 
metastable well, the quantum Kramers rate. The proposed formula works well in the current- voltage 
characteristic of a single Josephson junction. 
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The electrical conduction is a typical non-equilibrium 
phenomenon. Two theories are known. One is the Green- 
Kubo-Nakano linear response theory JjJ and the other 
is the Landauer theory 0. The former has a general- 
ity but often technically challenging. The latter, on the 
other hand, lacks a generality and limited to mesoscopic 
systems so far such as quantum wires, but is intuitive, 
simple and practically quite useful. 

The system Landauer considered is composed of a sam- 
ple which is connected by two leads to the electron reser- 
vors with chemical potentials /ii,2- The voltage applied is 
V = —fJ,2)/e. The electrons which contribute to the 
net current are those electrons with energies between fii 
and fi2- Write the energy of electrons with wave number 
k as E(k), then it contributes —e(dE/Tidk)T(k) where 
T{k) denotes the transmission rate. Assuming that any 
electron supply and removal from resevors are immedi- 
ate, the net current is simply, including factor 2 for up 
and down spins, 
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where fii ^ — E(ki^)- When fi\ r 
const = T, one gets the conductance 
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The purpose of the present paper is to propose a 
Landauer-like formula for calculating a particle's escape 
rate out of a dissipative metastable well, the so-called 
quantum Kramers rate. This well-defined and universal 
problem is ubiquitous in quite a wide range of physics 
and attracted much attention over decades Hf). The 
latest review article discussed a seemingly ultimate is- 
sue of full dynamical analysis in the presence of time- 
dependent external driving force || . In spite of all these 
quite intensive studies, however, the issue appears to be 
still open. The situation looks somehow similar to the 
problem of electrical conductance. As mentioned above, 
we have the linear response theory which is exact but 
often technically challenging and the Landauer theory 
which is limited but intuitive, simple and has seen a re- 
markable success in mesoscopic systems. It may be worth 



pursuing an analogous simple formula for the quantum 
Kramers rate. 

First note that (|lj) is a sum of individual contributions 
(dE /hdk)T(k) which is nothing but a right going flux. A 
corresponding quantity in the escape rate is the proba- 
bility of finding particle moving rightward at the top of 
the potential barrier. Adding all such contributions, we 
might reach a formula 
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Tr [e-^0(p)] /Z = ^<S(x)\p\> 
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where Z is the partition function of the metastable well, 
J- = 5(x)p/M with p — —ihd x is the flux at the barrier 
top and 0{p) is the unit step function. The formula (||) 
is nothing but the quantum transient state (QTS) theory 
long known in chemistry for calculating molecular reac- 
tion rates || . Historically, the QTS theory first appeared 
as an approximation to a formally exact expression from 
scattering theory H , 



Vexact = Re < TV > 
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where Re means to take the real part and the operator 
V = lim e lHt / h Q{p)e- lHt/n (5) 

t — >oo 

projects onto all states that have positive momentum in 
the infinite future. The exact expression (Q) can be ma- 
nipurated to give H 



dt < F(0)F(t) > 
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where J-(t) is the Heisenberg picture of the operator T . 

The formula (^) takes precisely the form of the Green- 
Kubo-Nakano formula for the electrical conductance. 
The QTS theory was obtained from the exact one (Q) 
by simply replacing V — > &{p) thereby neglecting all the 
dynamical processes. The QTS formula (||), however, is 
not well-defined quantum mechanically because the op- 
erator T is ambiguous. The necessary regularization of 
such operators as T is done by the Weyle rule |j| . Let us 
write (||) as 
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Now the Weyle rule gives a quantum operator from a 
classical one, J- C i(p,q) — > J- op (p,q), in coordinate matrix 
representation H, 



<C x | x ^> 



Applied to the last term in ([?]), we have 
< x'\S(x)\ P \\x > = 1 y dpe-*^*-*'^"!?^ (^4^ 

Putting (||) into (0) and integrating in parts gives 
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v ^ = irhpF- 1 d <--) d - < X \^ H \ ~ x > ( 10 ) 
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Some comments are due on the Weyle rule. First, it 
is not exact. For example, if one proceeds with this rule 
for the propegator e ~ lHt/h with H = p 2 /2M + V(x), the 
Weyle rule gives only a short time approximation . In 
general, as one can easily see from (^), the more compli- 
cated the functional forms of coordinate and momentum, 
unless they are functions of coordinate alone or momen- 
tum alone, the more dubious results are obtained by the 
Weyle rule. Back to the propagator, because it is already 
Hermite, the classical and quantum expression take the 
same form. Thus the failure of the Weyle rule resides 
in an inability of doing proper commutation calculations 
for complicated mixtures of coordinate and momentum. 
In fact when the starting classical form is simple enough, 
e.g. xp, then the Weyle rule is known to give a correct 
answer, e.g. (xp + px)/2. It is certainly embarrasing 
that we still do not know how to construct correct quan- 
tum mechanical operators from classical quantities when 
they are complicated enough. Second, a connection to 
the more familiar Wigner distribution function should 
be mentioned. In fact we have for any operator O, 

<0>= / dxdx' < x\e~ 0H \x' >< x'\O op \x > (11) 



Using the Weyle rule (|8|) for the last term in ( |11[ ) and 
changing the variables as x — x' — r, 
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q, we have 

0>= J f dpdqO cl (p,q)W(p,q) (12) 



where the Wigner distribution function W(j>, q) is the 
Wigner representation of the density matrix || , 

W(p,q) = ± f dre-^' n <q+ r -\e- m \q- r -> (13) 



The expression (|12|) should not be confused with the for- 
mally exact expression 
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where Ow is the Wigner representation of the operator 
O, 



W = j dre' tpr / h <q+ r -\O op \q 



> 
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It is a simple exercise to check that Ow reduces to O c i if 
the term in (|l5|) is approximated by the Weyle rule (^). 
Third, as is seen in @ and (|To|), the term 8(q) — 5( x ^-) 
results in an antiperiodic paths < x\e~@ H \ — x > in 
the expression (|io|). The nonlocality, an integral over x 
rather than just a contribution at the barrier top x = 0, 
however, arises from the step function O(p). It is noted 
that some nonlocal contribution is certainly negative re- 
flecting the fact that the Wigner distribution function 
is a quasi-probability function. It is also noted that 
the nonlocality does not show up when the QTS for- 
mula is treated by a semiclassical approximation (see e.g. 
Eq(2.22) in the second paper in M). 

We have critically reviewed the QTS theory ( |To| ) as a 
candidate for the Landauer-likc formula for the quantum 
Kramers rate. Two major approximations involved are 
(A) static or transient-state approximation and (B) the 
Weyle rule (||). The step function 0(p) or equivalently 
\p\ under the trace operation is found to be a dangerous 
operator for the Weyle procedure, leading to nonlocal- 
ity and associated possible negativity. Our idea in this 
paper is to get rid of this dangerous operator by taking 
square and square root. We thus propose a Landauer-like 
formula 



vli = 7^V< S(x)p 2 > Jl 
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where L is the size of the potential well. The factor L 
and the square root operation arises dimensionally and 
the fact that T is a density operator containing S(x). The 
proposed formula also makes the Kramers rate calculable 
from the local, at the barrier top, particle state. The fac- 
tor 1/2 takes into account only half of the contribution, 
outgoing one. It may be worth emphasizing that the p 2 
operator rather than 0(p) or \p\ would make the Weyle 
rule less dubious. Now, proceeding as before, we have 
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which should be compared with the QTS formula (|Io|). 

Note that the proposed formula (jl7|), unlike the QTS 
formula pc|), is manifest local. In fact the difference be- 
tween the QTS formula and the proposed one resides in 
a quantum fluctuation. To see this, consider 



< 8(x)\p\ > 2 = ^2 < n\e 

n,n f 



fm 5(x)\p\ 
n >< n'\e~ l3H 5(x) 



\p\\n>/Z 



2 



where the plane wave complete set may be used for 
{n}. Neglecting a quantum fluctuation, the term [...] — > 
<5« n' < n\S(x)\n' >, we reach 



< 5(x)\p\ > 2 



< S(x)p 2 > /L 
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We have tested the proposed formula ( Jig ) versus the 
QTS formula (|) for the I — V characteristic of the sin- 
gle Josephson junction (JJ). The resistively shunted and 
current biased single J J is described by the action for the 
superconducting phase difference between the two super- 
conductors lol. 
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where C is the capacitance, Rt the shunt resistance, I 
the bias current and Ej the Josephson energy. Thermal 
over-the-hill motion or quantum tunneling of the phase 
gives rise to a voltage V = Ucb/2e (the Josephson rela- 
tion) with = [cf),H]/ih = (2e) 2 n/hC where n = —id^ 
is the number of the Cooper pairs. So the flux operator 
here is 



T = 



(2e) 2 
hC 
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where 0o denotes a maximum point of the potential U(<p). 
Taking into account the backward flux by a detailed bal- 
ance, we have 
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After some manipulations, one can express the average 
in (21) as 



(5(0 - o )n 2 > = Z/Z Q 
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where 0(0) = and M§K) = <j>' ■ 

The path-integral (^) with the action S given by jl9| ) 
can be evaluated precisely by the cluster transfer ma- 
trix (CTM) method E2). In the present problem, how- 
ever, the dimensionless junction conductance g = R q / Rt 



where Rq = h/4e 2 — 6.45 Kil is the quantum resistance, 
is at g < 1, and from the study of the single electron box 
which has a similar action as ([H]), this regime of g can 
be accurately handled by the 1-cluster TM method 

For a fixed g, the temperature dependence of the re- 
sistance changes from insulator-like behavior ^ < to 

superconductor-like > with incresing ratio Ej / Ec 

2 

where Ec = Repeating the calculations for differ- 
ent g, we can thus map a superconductor-insulator (SI) 
phase diagram in the Ej / Ec — g plane. Fig. 1 shows a SI 
phase diagram at T=80 mK. The corresponding exper- 
imental results are denoted by open circles (S-like) and 
solid circles (I-like) |TJJl^]. Our result is the open di- 
amond, the band-theory result the thick solid line [fl6|| , 
and the QTS result is denoted by triangles. While there 
is an issue of temperature dependence concerning the 
phase diagram |l5} ], the agreement between the previ- 
ous theory and experiments with the proposed method, 
not with the QTS theory, will be evident. A major dis- 
agreement between theory and experiment is for some 
data points near g = 2.8 fig ]. However a similar phase 
diagram experimentally found for the 2D JJ arrays with 
similar parameter ranges for g, Ej / Ec and T is bounded, 
Ej/E c < 0.5, and g < 0.5 (cf. Fig 3 in 0). The above 
data near g = 2.8 is currently mysterious. 

In conclusion, we have proposed a new formula for 
calculating the quantum Kramers rate which may cor- 
respond to the Landauer formula for electrical conduc- 
tance. The single Josephson junction for which the pro- 
posed formula was tested is quite general: It contains 
quantum fluctuation, dissipation and external bias. The 
reasonable outcome of the proposed formula in compari- 
son with experiments and the previous theory, although 
the latter is not free from criticism [fl3"[ , may be encourag- 
ing a further testing of the proposed method in a variety 
of physical systems. On the other hand, we must also 
admit that the new formula was derived in an attempt 
for technical improvement over the QTS theory. A con- 
templation of deep physical argument leading directly to 
the proposed formula is much desired. 
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FIG. 1. Phase diagram of the single Josephson junc- 
tion. The phase boundary lies between the insulator- 
like (solid circles) and superconductor-like (open circles) 
samples experimentally found in [14,15]. The thick line is 
the band theory [16] at T = 80 mK. The triangle is the 
QTS theory and diamond is due to the present theory at 
T = 80 mK. 
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